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ABSTRACT 

We investigate the thermal response of the atmosphere of a solar-type star 
to an electron beam injected from a hot Jupiter by performing a 1-dimensional 
magnetohydrodynamic numerical experiment with non-linear wave dissipation, 
radiative cooling, and thermal conduction. In our experiment, the stellar at- 
mosphere is non-rotating and is modelled as a 1-D open flux tube expanding 
super-radially from the stellar photosphere to the planet. An electron beam is 
assumed to be generated from the reconnection site of the planet's magneto- 
sphere. The effects of the electron beam are then implemented in our simulation 
as dissipation of the beam momentum and energy at the base of the corona where 
the Coulomb collisions become effective. When the sufficient energy is supplied 
by the electron beam, a warm region forms in the chromosphere. This warm 
region greatly enhances the radiative fluxes corresponding to the temperature of 
the chromosphere and transition region. The warm region can also intermittently 
contribute to the radiative flux associated with the coronal temperature due to 
the thermal instability. However, owing to the small area of the heating spot, the 
total luminosity of the beam-induced chromospheric radiation is several orders 
of magnitude smaller than the observed Ca II emissions from HD 179949. 



Subject headings: MHD — methods: numerical — stars: atmospheres — plane- 
tary systems 
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Introduction 



Hot Jupiters are Jupiter-mass planets located within ~ 0.1 AU or less from their par- 
ent stars. Because of the close proximity to the parent stars, hot Jupi ters have been ex- 
pected to be able to influen ce their stellar companior is via magnetic ( 



Rubenstein fc Schaeferll2000l ) and/or tidal interactions (ILin et al.lll996l : iJackson et al 



Cuntz et al 



2000 



2009 



Pfahl et al. 20081 ). The observations of Ca II H & K lines from a number of stars harbor- 
ing a hot Jupiter have suggested that the chromospheric activities, characterized either by 
l ine intensity or shor t -time variab ility, sometimes correlate with the orbits of their planets 
( IShkolnik et al.ll2003l . |2005| . |2008| ). These phenomena can be modelled as a hot spot or a 
more "variable" region, despite residing in the chromosphere, following the planet's orbital 
motion with a phase difference. In particular, the observations carried out in 2001, 2002, 
and 2005 imply that a hot spot on HD 179949 persistently leads the planet by ~ 70° with 
the intensity of ~ 10^'' erg/s in Ca II emissions. The similar phase lead of a variable region 
in optical has been suggested by the MOST satellite photometry for the hot- Jupiter host 
star r Bootis (IWalker et al.ll2008l ). Since the planet-induced stellar activities occur only once 
during one orbital period, the origins of this "spot" have been attributed to magnetic rather 
than tidal interactions. 

One of the commonly adopted scenarios to describe the star-planet magnetic interactions 
is the magnetic interactions between Jupiter and its Galilean satellites (see Zarka 2007 for 
a review). In this scenario, the orbital motion of the Galilean satellites relative to Jupiter's 
magnetosphere taps the orbital energy of the satellites at a rate that depends on detailed 
modelling on the magnetic interactions. The interactions can be classified into two types: the 
unipolar interaction with an unmagnetized satellite such as lo and the magnetic reconnection 
with a magnetized satellite such as Ganymede. The energy is then transported by the Alfven 
waves and/or by a fast electron beam along the field lines from the satellites to Jupiter's 
surface where the energy is dissipated, thereby explaining the satellite-induced emissions 
from Jupiter. In the case of star-planet interactions, this picture has been modified to take 
into consideration stellar winds along open field lines or to allow for a large static magnetic 
loop connecting the star and the planet. As a result, the phase differences between the 
stellar "spot" and the planet are expl ained by the time la g due to the Alfven travel time 
to the star in the Alfven- wave model (jPreusse et al.ll2006l ). or by the larg e magnetic loop 
having a geometry across longitudes of the star in the electron-beam case Jlanzalboosh . A 
three-dimensional resistive magnetohydrodynamic simulation was performed to study how 
the magnetic field- aligned current can develop from a hot Jupiter (jPreusse et al.l 120071 ). 



However, how the stellar atmosphere thermally responds to any energy injection from 
the planet so as to generate the chromospheric emissions of > 10^^ erg/s remains elusive. The 
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magnetic energy at the magnetopause of a hot Jupiter has been estimated t o be insufficient 



to supply the energy rate of > 10^^ erg/s fIShkolnik et a. 
of the stellar surface field is ~ a few Gauss (ICatala et al. 



20051 : IZarkall2007n if the strength 



20071 ). Furthermore, if this energy 



dissipation rate arises entirely from the orbital energy of a hot Jupiter, most hot Jupiters 
would plunge into their cent ral stars in a few billi on years, a timescale comparable to the age 
of these planetary systems (jShkolnik et al.ll2005l ). However, the gas in a stellar atmosphere 
is certainly not quiescent but is fluctuating with the free energy t hat may be l i berat ed as 
a hot Jupiter enc ounters the turbulent stellar fields along its orbit. ICuntz et al.l (120001 ) and 
Saar et al.l (120041 ) took into account the energy contribution from the stellar macroturbulence 
velocity and compared the stren gth of the planet -star interactions relative to each other for a 
number of hot- Jupiter systems. iGu et al.l (120051 ) postulated that most of the planet induced 
emissions may result from stellar turbulent energy to reconcile the energetic problem. To 
further examine this possibility, a more realistic stellar atmosphere model implemented with 
an energy equation involving stellar turbulence is required. 

In the case of the solar atmosphere, the coronal fields that can reach the typical or- 
bits of hot Jupiters are the open magnetic fields emanating from the coronal holes. The 
dynamical features of the open fields are the non-linear fluctuati ons in the corona a r id the 
high-speed winds (800 km/s at ~ 1 AU). To explain these features. ISuzuki fc Inutsukal (120051 ) 
(hereafter SI05) introduced a 1-dimensional magnetohydrodynamic simulation with radia- 
tive cooling and thermal conduction. In this simulation, the authors have self-consistently 
treated the transfer of mass/momentum/energy by solving the magnetic waves propagat- 
ing from the photosphere to the interplanetary space. The open field lines are a ssumed to 
expand super-radially from the photosphere to the corona (IKopp fc Holzerlll976l ). which is 
consistent wi th the spectropolarim etric measurements of magnetic structures in the Sun's 
polar region ( Tsuneta et al.ll2008l ). The heating and acceleration of the gas in the coronal 
holes are achieved by the dissipation and the pressure of the non-linear magnetic waves. The 
model can explain the structure of the solar atmosphere (i.e. consisting of the chromosphere, 
the thin transition region, and the corona) as well as the high-speed winds from the coronal 
holes. 

In this paper, we employ the magnetohydrodynamic simulation described above and 
model the open fields as a 1-dimensional flux tube from the stellar photosphere to the hot 
Jupiter. The aim is to conduct a numerical experiment as a starting point to study the 
thermal properties of the stellar chromosphere, transition region, and corona in response to 
an energy injection. Our study concerns HD 179949, but it is also guided by information 
about the solar wind due to the fact that various properties of HD 179949, especially in 
regard to its chromosphere, corona and wind, are poorly constrained. In §2, we describe 
the stellar atmosphere model in the presence of an inward-propagating electron beam. The 
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numerical results are presented in §3. The paper concludes with a summary and discussion 
in 54. 



2. Model Description 

We consider a magnetized hot Jupiter revolving in a circular orbit on the equatorial 
plane of its parent star. As the planet orbits through the open flux tubes emanating from the 
stellar surface, magnetic reconnections occur and generate plasma jets propagating inwards 
along the flux tubes to the parent star. In the following subsections, we shall describe the 
model of electron-beam injection and the stellar atmosphere associated with the flux tube. 



2.1. Electron-Beam Injection 



Since HD 179949 has been the canonical planetary system highlighted by most previous 
studies for magnetic interactions, we adopt the parameters of HD 179949 as an illustrative 
example for our main study. In other words, we focus on a hot Jupiter orbiting at the radial 
distance r = 7.8-R* around a central star of mass = 1.21Mq, radius i?* = 1.22i?0, and 
effective temperature Tes = 6168 K (IButler et al.ll2006l ). In terms of M^,, R^, and Tes, these 



stellar parameters are noticeably different from the solar values, although in the broad sense 
HD 179949 still constitutes a solar-type star. 

The typical field strength at the plane orbit is i? ~ 0.1 — 0.01 G, for the average radial 
field strength ~ 1 — 10 G at the stellar surface (i.e. B oc r^^). At r ^ IORq, the planet is 
located inside the Alfven radius of the star (see the next subsection) and therefore no fast 
MHD shocks form as the stellar winds encounter the planet's magnetosphere (cf. Zarka 2001; 
Ip et al. 2004; Preussue et al. 2005). Assuming a dipole field for the planet's magnetic field, 
we can estimate the radius of the magnetopause Rmp at which the stellar and the planet's 
fields are balanced: 



Rmp = ^.U^Jtj 



AMR.J 

\Rj 



Bp 

Tg 



Tg 



IORq 



R* 

Rq 



1 1/3 



where Bp and 5* are the average surface fields of the planet and the parent star respectively, 
Rj is the radius of Jupiter, and a is the semi-major axis. As can be clearly seen from the 



above equation, R 



■mp 



<< a 



IORq. Therefore the magnetic field strength Bmp 



Bmp ~ 



at the 
0.1-O.OlG. 



magnetopause is almost equal to the stellar field strength at r = a; i.e. 

We consider a jet originating from the reconnection sites at the magnetopause. The 
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energy flux liberated by the reconnection is on the order of {B"^^/ An)vreii where Vrd is the 
speed of the planet relative to the stellar fields. For a slowly rotating solar-type star, v^ei is 
cl ose to the orb ital speed of the planet ^JGM^^^^ = 1.38 x 10^ cm/s (cf. Vrei ^ 300 km s~^ 



in IZarkal (120071 )). Thus, we arrive at a rough estimate of the liberated energy flux 



A fraction of this energy is converted to the kinetic energy of the reconnection jets. The speed 
of the jets is on the order of the local Alfven speed -Bmp/ ^/Amp ~ 900km s^^ (oMg) ( lo-^igcm-a 
where p = 10~^^ g/cm^ is the typical mass density at r ^ IORq (e.g. Suzuki & Inutsuka 
2006). We assume that the same energy per unit mass is converted to the kinetic energy 
of the electrons. Therefore the electrons are moving faster than the ions by a factor of the 
square root of the ion-electron mass ratio (~ 40). As a result, f t ~ 3 x 10^ km s~^ is a 
typical speed of the electron beam. 

We denote the momentum flux and the energy flux of the electron beam as Pb and Ff,, 
and their initial values as Pbi and F^i, respectively. If we neglect thermal fluctuations of the 
beam particles, the initial momentum flux and energy flux of the electron beam are 

Pbi = Pbivl, (3) 
Fbi = -Pbivt, (4) 

where pbi is the initial beam mass density. In the present numerical experiment, we restrict 
ourselves to the cases in which Fbi = 10^ and 10^ erg cm~^s"^. We adopt a constant 
f f, = 3 X lO'^km s~^, and vary p^ for different Fbi. Note that according to eq. ([2]) and 
the parameters of the HD 179949 system, Fbi = lO^erg cm~^s~^ corresponds to Bmp ~ 
0.005 — 0.01 G and hence 5* ~ 1 G, and the larger energy flux Fbi = lO^erg cm"^s~^ 
corresponds to Bmp ~ 0.05 — 0.1 G and thus the stronger stellar field B^ ^ 5 G. 

We assume that the electron beam initially streams freely along the flux tube at a 
constant Vb- As they approach the parent star, the beam particles start to interact with the 
dense background ions through Coulomb collisions and become more concentrated due to 
the flux-tube convergence (see the next subsection). The collisions lead to the heating of 
surrounding plasma. The heating becomes efficient when the mean free path of the beam is 
comparable to the local density scale height Hp. The mean free path of an electro n colliding 



with a pool of thermal ions is given by /mfp ~ 3 x lO'^km ( gxio^ms-i )^ ( io9cm-0 ^ (iBraginskii 



19651 ). where n = lO^cm"^ is the typical density at the upper transition region or the lower 
corona (e.g. see SI05, or refer to the density profile in Figure 1 to be discussed in §3). Roughly 
speaking, for Vb = 3 x 10^ km/s, /mfp < Hp when n > lO^cm"'^. Therefore, we anticipate 



-1/2 
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the inward-propagating electron beam to heat up ambient media from the lower corona to 
the upper chromosphere. This allows us to assume that the incoming electron beam starts 
dissipating from rmax = l-l-R© (lower corona) to r^^^ = l.OOli?0 (upper chromosphere). 

We model the energy flux of the beam to decrease inwardly according to 

F,^( "'""^'^ V. (5) 



where is a parameter that describes the spacial distribution of the heating, k = 1 corre- 
sponds to constant volumetric heating; namely, the heating rate per unit mass is higher in 
the upper region (i.e. corona) than in the lower region (i.e. chromosphere). If = 0.1, the 
heating rate per unit mass is more uniformly distributed. In this work, we adopt k = 0.1 
to resemble the situation of a constant beam-heating rate per unit mass. We also assume 
that the momentum flux of the beam Pf, dissipates in the same manner as that described by 
eq. ([5]) for the energy flux. This implies that although the beam velocity Vb does not decay 
in our dissipation model, the beam density pb declines, meaning that more and more beam 
electrons have been transformed into thermal electrons as the beam progresses downwards 
in the dissipation region. 

The beam heating at the footpoint of one open field line occurs when the planet's 
magnetosphere is crossing the field line. Hence, by means of eq. (Q, the beam heating 
proceeds on the timescale 



heam = ^ 80 mius ( 



\RjJ V138km/s 



-1 



1/3 



(6) 



2.2. Stellar Atmosphere Model 

The stellar atmosphere model is based on the 1-D magnetohydrodynamic simulation 
with radiative cooling and thermal conduction in an open flux tube (SI05). In the original 
simulation, the heating is given by the nonlinear dissipation of the Alfven waves excited by 
the granulations at the photosphere. In the case oi = 1 G, we set a rms average amplitude 
< dv±_ >~ 1.8 km/s at the photosphere, which is estimated from the scaling with the surface 
convective flux (Suzuki 2007) from the Sun (SI05). In the case of the stronger field i?* = 5 
G, the larger rms velocity fluctuation < dv^ >~ 3.6 km/s is used for the experiment, which 
is twice the value for the -B* = 1 G case. We neglect the effect of rotation for simplicity in 
order to focus mainly on the energetics influenced by the electron beam. 
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When an inward-propagating electron beam described in the above subsection is taken 
into account, the momentum and energy equations in a 1-D open flux tube are modified to 
(cf. SI05) 

and 

+^|:(-VFc) + ^i^r^fF.) + = 0, (8) 

where p, V, p, e, and B are the density, velocity, pressure, specific internal energy, and 
magnetic field strength, respectively; the subscripts r and ± denote radial and tangential 
components, respectively; d/dt and d/dt denote Lagrangian and Eulerian derivatives, re- 
spectively; G and are the gravitational constant and the stellar mass, respectively; Fc is 
thermal conductive flux; qr is the radiative cooling and / is a super-radial expansion factor. 
We assume that the super-radial expansion (in addition to the radial expansion oc r^) is 
a factor of 240 and 480 from the photosphere to 2Rq in the cases of = 1 and 5 G, 
respectively. It then follows that to give -B* = IG and 5G, the radial magnetic field strength 
at the footpoint of the flux tube at the photosphere Bph is 240G and 2400G, respectively. 

Note that on account of the convergence of the flux tube toward the central star, the 
energy flux of the incoming electron beam increases due to the areal focusing. In the case of 
B^ — IG, the flux tube converges by a factor of ~14,600 (a factor of ~ 61 from the radial 
convergence and another factor of 240 from the super-radial convergence) from 7.8-R* (the 
planet's orbit) to li?*. On the other hand, the converging factor of the flux tube in the 
B^ = 5G case is twice as large. In the absence of dissipation, the energy flux of the beam 
increases by the same factor as the converging factor of the flux tube. In contrast, the same 
enhancement effect of the electron beam pressure on stellar thermal ions is almost negligible 
in our calculations because of the small mass of an electron. 



3. Results of the Numerical Experiment 

We conduct the numerical experiment to study the effects on the stellar atmosphere 
due to the electron beam with the initial energy ffuxes for the two B^ cases as described in 
the preceding session. The experiment is flrst carried out in the absence of the beam until 
the simulated atmosphere attains a quasi-steady state. The beam heating is then added 
afterwards. 
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Figure 1 shows the snapshot structures of the stellar atmosphere for = 1 and 5 G 
denoted respectively by "small" and "large B and dv" , and compares the results with and 
without the electron-beam injection. In the B^ = IG case, the results for Ffyi = 10^ erg 
cm~^ s~^ (blue curves) are almost identical to those for the no-beam case (green curves). In 
contrast, the atmospheric structures change noticeably in the larger B^ (= 5G) and < dv± > 
(= 3.6 km/s) illustrated by the red and black curves. In the absence of the beam 

heating, the stronger dissipation arising from the larger < dv± > raises the chromospheric 
temperature and therefore causes the chromosphere to evaporate and expand (see the red 
curve in the middle panel). The local density scale height increases and the density drops 
more slowly with r. Therefore, the locations of the transition region lie at higher altitudes 
in the larger < dv± > case. A hot and dense region forms in the chromosphere, which we 
term as "a warm region" (i.e. T > lO^K, see the bottom panel). It then becomes difficult to 
further heat up the warm region to the c oronal temperature (i.e. > 10^ K) as a result of the 



efficient radiative cooling at T < 10 K (ILandini &: Monsignori-Fossilll990l ). However, when 



the beam energy of Fu = 10'^ erg cm~^ s~^ is added, the chromosphere is further heated up 
and thus an even larger warm region can form at the location of r = 1.001 — 1.01-R^, (see the 
black curve in the bottom panel). The black curve in the middle panel indicates that the 
density is on the average much higher in the beam-heated region. Consequently, the heating 
rate per unit mass becomes smaller and the temperature of the coronal region at r > 1.01i?^= 
then becomes lower in the case of Fbi = 10^ erg cm~^ s~^ than in the case of Fbi = (see 
the bottom panel). The sharp transition region disappears. 

The top panel of Figure 1 shows the radial profile of the stellar wind velocity Vr- The 
winds in the outer region (r > 2_R*) are faster in the large 5* case. 5* equals Bph/ f, which 
in fact specifies the fiux tube properties. The larger the factor Bph/ f is, the r nore the wave 



energy dissipates in the outer region of the atmosphere, leading to faster winds (IKojima et al. 



20051 : ISuzukil l2006l ). However, the winds on the average are not significantly amplified by 
the beam heating, as indicated by the large overlap between the red and black curves in the 
plot. The wind speed at the planet's orbit is ~ 300 km s~^, which is lower than the Alfven 
speeds ~ 500 and 1000 km s~^ there in our cases for B^, = 1 and 5 G. That is, the planet 
lies inside the Alfven radius. 

Having described the snapshot structures, we should note that the thermal properties of 
the stellar atmosphere actually fiuctuate with time owing to wave propagation and dissipa- 
tion. As a result, the location and thermal properties of the warm region fiuctuate with time 
as well. Figure 2 shows the evolution of the radiative fiuxes arising from different temperature 
components in the beam and no-beam cases. In general, all of the radiative fluxes fluctuate 
with time. Nevertheless for the case of 5* = IG and < v± >= 1.8 km/s, the electron beam 
of Fbi = 10^ erg cm~^ s~^ gives rise to only tiny effects on the radiations compared to the 
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no-beam expected from the snapshot structures shown in Figure 1. In contrast, the 

case of F^i = 10^ erg cm~^ s^^ corresponding to -B* = 5G exhibits noticeable differences from 
the no-beam case, which can also be expected from the snapshot structures in Figure 1. The 
radiative fluxes from the hot chromosphere (7000-20000K) and the transition region (20000- 
5 X lO^K) are greatly enhanced by the warm region in the vicinity of r = 1.001 — l.Oli?*. 
Since the warm region is on the average denser than the usual chromosphere and transition 
region in r > 1.003-R,. (see Figure 1), the radiative flux corresponding to the chromospheric 
temperature is on average increased by a factor of ~ 10 and the radiative flux associated 
with the transition-region temperature is intensified by a factor of ~ 100. 

The bottom panel of Figure 2 shows that the beam of Fhi = 10^ erg cm~^ s^^ associated 
with = 5G is able to enhance the radiation from the hot gas of T > 5 x lO^K occasionally 
by a factor of > 100. This radiation is mainly from the warm region that intermittently 
develops around r = 1.001 — 1.003-R*. With this beam heating, the temperature of the 
warm region can sometimes go up and do wn between lO^K and lO^K. This temperatur e 



variability is due to the thermal instability (ILandini fc Monsignori-Fossilll990l : ISuzukill2007l ). 



The electron beam is energetic enough to continuously heat up the warm region, while the 
wave dissipation heats it up in a stochastic manner. A small change of the wave heating rate 
triggers violent fluctuations of temperature in the thermally unstable regime of 10^ < T < 
lO^K. 

As is shown in Figure [H the thermal response of the stellar atmosphere to the onset of 
the beam-heating to evolve to the hot state is nearly instantaneously. In §2.11 we estimated 
the duration of the electron beam as ~ 100 min, which is rather short in comparison with 
the simulation time presented in Figure [2l To estimate the area of the hot spot, it is equally 
important to examine how long the hot atmosphere can be sustained after the electron beam 
passes by. Figures [3] & H] show the results for the -B* = 5G case but the beam is switched 
off at t = 109 min. Figure [31 which presents the snapshot structure at 219 min, namely 
110 min after the beam stops, illustrates that while the warm region in the chromosphere 
is still present as shown by the temperature profile, its density has dropped to the original 
lower-density level. Figure H] shows that the radiative fluxes from the hot chromosphere (top) 
and transition regions (middle) have declined for ^ 50-60 minutes since the beam is switched 
off, which indicates that the atmosphere takes only a fraction of the heating time to revert 
to the normal state in the absence of the electron beam. Therefore, the size of the hot spot 
estimated in terms of the projected area of the planet's magnetosphere along a flux tube 
onto the chromosphere is a reasonable approximation. 



Shkolnik et al.l (120051 . l2008f ) found that in the case of HD 179949, the planet induced 
Ca II flux averaged over the stellar disk is ~ 1.5 x 10^ ergs cm~^ s~^, which amounts to 
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an increase of a factor of about 1.04 compared to the non-planet induced Ca II emissions. 
Hence the Ca II emissions with no planet induced component is ~ 1.5/0.04 x 10^ = 3.7 x 10^ 
erg cm~^ s~^, which is on the similar order of the value given from the model shown by 
the red-dotted line in Figure O In this sense, the model for 5* = 5G and < dv\_ >= 3.6 
km/s may mimic a chromospheric condition similar to that of HD 179949, although our 
cooling function (/r for th e chromospheric radiation is based on the observation of the Sun 
(lAnderson fc Athaylll989l ). Having obtained the radiative flux ~ 10^~^ erg cm~^ s~^ from 
the hot chromosphere in the case for 5* = 5G and < dv± >= 3.6 km/s, we can estimate 
the luminosity of the hot spot. In our model, the cross-sectional area of the converging flux 
tube in the chromosphere is about 8000 times smaller than the planet's magnetosphere. In 
other words, the area of the hot spot in the chromosphere is given by 

A.,romo ^ 10''-'' ( cm\ (9) 



5Rj 

which when multiplied by the chromospheric flux gives the luminosity of the chromospheric 
hot spot ~ xo^^~^^ erg/s. This total chromospheric emission is still 2-4 orders of magnitude 
weaker than the observational Ca II emissions for HD 179949. 



4. Summary and Discussions 

By conducting a numerical experiment, we study the thermal response of the atmosphere 
of a solar-type star to the dissipation of an injected electron beam at the coronal base. The 
experiment is carried out based on the framework of the 1-D magnetohydrodynamic simu- 
lation by SI05 with non-linear wave dissipation, radiative cooling, and thermal conduction. 
We assume that the magnetic stress due to the orbital motion of the planet relative to the 
stellar coronal flelds generates an electron beam, which in turn funnels along the stellar open 
fleld lines to the central star. As the beam travels inwards, the energy flux of the incoming 
electron beam is intensifled by the areal focusing of the super-radially converging open flux 
tube. 

We use the stellar parameters of HD 179949 as an illustrative example but ignore possible 
magnetic properties arising from its stellar rotation. When the average stellar fleld at the 
photosphere is about 1 G and the average amplitude of the wave velocity < v± > is about 
1.8 km/s, the stellar atmosphere is not considerably altered after the beam dissipation is 
turned on. In contrast, when =5 G and < v± >= 3.6 km/s, we flnd that a warm 
region forms in the chromosphere. The warm region becomes substantially hotter and denser 
once the electron-beam heating is switched on. As a result, the beam-intensifled warm 
region enhances the chromospheric radiative flux by a factor of ~ 10, and the radiative flux 
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corresponding to the temperature of the transition region by a factor of ~ 100. The warm 
region can also intermittently contribute to the radiative flux associated with the coronal 
temperature by a factor of ~ 100 due to the thermal instability (iLandini &: Monsignori-Fossi 
1990l ). In other words, the planet-induced radiations are not perturbations in the local region 
of the hot spot compared to the normal state of the stellar emissions. However, owing to 
the small area of the heating spot, the total luminosity of the beam- induced chromospheric 
radiation is 2-4 orders of magnitude smaller than the observed Ca II emissions from HD 
179949. 

The energetics of the planet-induced emissions becomes a more serious problem in our 
numerical experiments when explaining the statistical results of X-rays from planet-host 
stars: stars with close-in giant planets are o n average more X-ra.y active by a factor ~ 4 
than those with planets that are more distant (IKashyap et al.ll2008l ). Since the typical X-ray 
luminosity from a solar-type stars is > ]^q27-28 erg/s, the planet-induced X-ray inferred from 
the statistical analysis is actually even stronger than the planet-induced Ca II emissions from 
HD 179949. Our simulation results show that an ~ 100 times enhancement in X-ray due 
occasionally to the thermal instability of the small warm region contributes only an even 
more negligible perturbation to the total X-ray emission, rather than being comparable to 
it. 

We note that our open-field model for the chromospheric emissions is different from 
those occurring on the Sun where the emissions come primarily from the solar plage regions 
in closed magnetic loops. Needless to say, our 1-D numerical experiment restricts ourselves 
to exclude any mechanical and thermal influence of the heating area on neighboring open 
fields and closed magnetic loops. As such, our results leave an open question as to whether 
the thermal instability or any other magnetic instabilities (e.g. Lanza 2008; Ishikawa et al. 
2008) can be further triggered around the beam-heated spot to liberate more energy. 

In our numerical experiment -B* has been taken to be 1 and 5 G, which is consistent with 
the field measurern ents via spectrapolarimetry for the similar spectral type dwarf r Bootis 
( ICatala et al.ll2007l ). However, the field strength inferred from the Stokes V observation may 
be underestimated due to the cancellation of circular polarization arising from the opposite 
directions of B^, along the line of sight. It is normally expected that a faster rotator and 
therefore a s tronger X-ray emitter may possess stronger B^, than the Sun (IRiiedi et al.lll997l : 
Giidell 120071 ). For instance, B^ of HD 179949 has been assum ed to be ^ 8-9 x so lar value 



by scaling with the X-ray flux (jSaar et al.ll2004l ). Furthermore, iMontesinos fc Jordan! (119931 ) 
related B^ and the filling factor fg to the Rossby number Ro (= Prot/Tc, where P^ot is the 
rotation period and Tc is the convection turnover ti me) of a dwarf star. In the case of HD 



179949, Prot has been suggested to be about 7 days (IWolf &: Harmanedl2004l : IShkolnik et al. 
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20081 ) and Tc ~ 4.287 days may be inferred from its color B— V=0.503 (Noyes et al. 1984). 
Thus following iMontesinos &: Jordan! (Il993l ). we obtain the average field Bfg = 37.6 G and 
the filling factor fg = 0.02. If we assume that these magnetic properties are contributed 
mainly from the open-field region of our model, then Bfg = 37.6 G is equivalent to 5* 
and the filling factor fg = 0.02 may correspond to the super-radial expansion factor 1/ f 
i n our experiment. In realit y, some of the contribution may come from closed- field regions 
(IMontesinos fc Jordanlll993l ). which introduces additional complexity. 



The other free parameter that governs our numerical results is the power spectrum of 
MHD waves. The stellar macrotu rbulent velocity Vmnr- HP 179949 is expected to be 2 



times larger than that for the Sun (jSaar fc Ostenlll997l : ISaar et al.ll2004l ). In our experiment 



we adopt a larger < dv± > for the stronger 5^, case. However, the correlation between Vmac 
and < dv± > is not tested in this work and how these velocities are associated with 5* 
is not modelled in our numerical experiment. After all, in view of all of the uncertainties 
and complexities mentioned above, our results serve only as the fiducial examples for future 
studies. The numerical experiment covering a broader parameter space coupled with more 
magnetic-field measurements will be essential to further diagnose the problem. 

The effect of the centrifugal force is not implemented in the present numerical experi- 
ment. While it is a reasonable approximation for thermally driven winds from a slow rotator 
like the Sun, the magneto-centrifugal winds play an equally ir nportant role in accelerating 
stellar winds for a star rotating ~ 10 tiraes fas ter than the Sun (IBelcher fc MacGregorlll976 



Washimi fc Shibatalll993l : iPreusse et al.ll2006l ). The centrifugal force of a 7 day-period star 
is ~ 1/5 that of a 3 day-period star, as the force is proportional to the square of the rotation 
frequency. The effect of the centrifugal force on the stellar atmosphere will be investigated 
in a future work. 

While most of the previous magnetohydrodynamic simulations and studies have focused 
on the planet's side for the magnetic interactions (Ip et al. 2004; Preusse et al. 2006, 2007; 
cf. Laine et al. 2008 for a young hot Jupiter), our numerical experiment makes an attempt 
to investigate how a stellar atmosphere down to the photosphere in the open field region 
responds to the dissipation likely from a hot Jupiter. The current model is simplified in such 
a way that we prescribe the energy dissipation and ignore stellar rotation. Nevertheless, the 
dissipation contributed by the planet is described only by the energy fiux Fb at the coronal 
base, meaning that Fh is not necessarily specific only to an incoming electron beam but 
can be in the form of other dissipative sources (e.g. damping of Alfven waves) with proper 
modification. Furthermore, our simulation lays the framework to extend the calculations for 
other types of dwarf stars if the corresponding magnetic wave amplitudes and power spectra 
are specified. 
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Fig. 1. — Snapshot structures of the radial velocity, density, and temperature for the different 
beam energy fluxes at 142 minutes after the electron beam is switched on. The green dashed 
and blue dash-dotted curves are the results of F^j = and lO^erg cm ~^s"^ for i?^, = 1 G, 
whereas the red dotted and black solid curves are the results of F^i = and 10"^ erg cm~^s~^ 
for = 5 G. Note that the vertical scale of the radial velociy is different from that of the 
density and temperature. 
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Fig. 2. — Comparison of the radiative flux arising from different temperature components. 
The top, middle, and bottom panels show the radiative flux from the gas at the temperatures 
corresponding to the hot chromosphere, transition region, and corona. The cases for B^, = 1 
and 5 G in the absense and presence of the beam heating are plotted with the same color 
coding as those shown in Figure [H 
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Fig. 3. — Comparison of snapshot structures of the density and the temperature at t = 219 
minutes for the case with the continuous beam of the energy flux = lO^erg cm^^ s^^ 
(black sohd curve) and the case with the beam switched off at 109 minutes (red dotted 
curve). 



- 19 - 




Beam stops 
Continuous beann 

_l 1 1 1 1 I L. 



_I I I I L_ 



10' 



10^ k 



T 1 1 1 1 1 1 1 — ITT 1 1 1 1 1 1 1 1 1 r 



5x1 o^k-' 




I:- 



50 



100 
time(min.) 



150 



200 



Fig. 4. — Comparison of the radiative flux for the case with the continuous beam (sohd) and 
the case with the beam switched off at 109 minutes (red dotted). The three panels show 
the radiation from the gas at the temperatures corresponding to hot chromosphere (top), 
transition region (middle), and corona (bottom). 



